Color-dependent conductance of graphene with adatoms 
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We study ballistic transport properties of graphene with a low concentration of vacancies or 
adatoms. The conductance of graphene doped to the Dirac point is found to depend on the relative 
distribution of impurities among different sites of the honeycomb lattice labeled in general by six 
colors. The conductivity is shown to be sensitive to the crystal orientation if adatom sites have a 
preferred color. Our theory is confirmed by numerical simulations using recursive Green's functions 
with no adjustable parameters. 

PACS numbers: 73.63.-b, 73.22.-f 



a 

E 

c 

o 
o 



> 
00 

o 



X 



Controllable deposition of adatoms known as chemical 
functionalization is an efficient way to alter graphene's 
electronic properties. Adatoms can change the orbital 
state of the functionalized carbon atoms from sp 2 to 
sp 3 configuration removing electrons from the conduc- 
tion band and transforming graphene into a semicon- 
ductor [HQ. For instance, fully hydrogenated graphene 
(graphane) is predicted to have a wide band gap of 3 — 14 
eV [3j in sharp contrast to clean graphene which has a 
gapless excitation spectrum 

A small concentration of strongly bound adatoms or 
molecules, such as hydrogen or CH3, is naturally present 
even in the cleanest graphene samples produced by ex- 
foliation. The effect of vacancies, which can be thought 
of as infinitely strong on-site potential impurities, is es- 
sentially analogous to that of adatoms. Such impurities 
support electronic states at the Dirac point, yielding res- 
onant scattering of Dirac quasiparticles |5l4l2j . 

From a theoretical point of view, the strong on-site 
impurities or vacancies preserve both time-reversal and 
chiral symmetries of the graphene Hamiltonian. The in- 
terplay of the symmetries prevents quantum localization 
at the Dirac point [5j and gives rise to a quantum criti- 
cal regime of charge transport at sufficiently high impu- 
rity concentration. The conductivity of graphene in this 
regime is predicted to take on a constant value, which 
depends on the distribution of adatoms among different 
sublattices of the graphene crystal [ll|. I n this Letter, 
we discover an even more subtle sensitivity of transport 
quantities, not only to sublattice but also to the type of 
impurity sites therein (described by the site's "color"), 
that can be clearly seen even in the presence of two scat- 
terers. 

The theory of quantum transport in disordered 
graphene is highly non-trivial in the vicinity of the Dirac 
point due to the breakdown of the quasiclassical approx- 
imation. An essentially exact approach, which is referred 
to below as an unfolded scattering formalism, has been 




FIG. 1: (Color online) Color scheme for vacancies or resonant 
adatoms. The impurity site is characterized by the phase 
0± = ±q + 2-7tc/3, where ± refers to the sublattice (A or 
B) and c = —1,0,1 denotes the colors (red, green, blue). 
The color scheme depends on the transport direction x. The 
phases at A and B sites connected by the vector So are equal 
for a — 0, while they differ most for a — n/6. 



proposed to tackle the problem in Ref. In this pa- 

per, we perform for the first time a detailed comparison 
of this formalism with numerical simulations based on 
the well-established recursive Green's function technique. 
In particular, we calculate both numerically and analyt- 
ically the conductance of graphene at the Dirac point in 
the presence of two vacancies. This quantity is shown 
to have a remarkable sensitivity to the "color" of the va- 
cated sites, which is determined by the phase of the Bloch 
wave-function and by the orientation of the crystal with 
respect to the transport direction. We further use the 
analytic results in order to construct the first terms of 
the color-dependent virial expansion of the conductivity 
with respect to vacancy concentration. 
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An essential in gred ient of the analytical approach de- 
veloped in Refs. [lOL |ll| is a single-impurity T- matrix, 
which is known for a variety of impurity types in 
graphene [13[ and carbon nanotubes [14j. Starting from 
the standard tight-binding model on the honeycomb lat- 
tice, one can relate the onsite potential Vi to the cor- 
responding T-matrix by the geometric series Ti = Vi + 
ViguVi + ViguViguVi + . . . , where g^ is the free particle 
Green's function taken between the lattice sites i and j. 
The key simplification occurs at the Dirac point (e = 0) 
since the Green's function element gu oc e In e vanishes 
in the limit e — > 0. The vanishing amplitude of return 
yields at e = a simple form of the T-matrix operator: 
Ti = Vi. Expressing this T-matrix in the valley-sublattice 
space of the Dirac Hamiltonian reveals a color scheme for 
the graphene lattice, which we shall describe below in de- 
tail. 

To be more specific we write down the tight-binding 
model of ideal graphene as 



Cr<5 s ) = e*(r), 



(1) 



where the sign factor £ r = ± specifies the sublattice of 
the r site (A or B),t fx 2.7eV is the hopping integral, and 
the index s takes on three values: 0, 1, and 2. We regard 
x as transport direction and orient the graphene lattice at 
an arbitrary angle a as shown in Fig. [TJ The sites of the 
honeycomb lattice are then connected by three vectors 
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s = 0,l,2, 
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where a is the distance between the neighboring carbon 
atoms. The two non-equivalent Dirac-point vectors in 
the reciprocal space, K = (47r/3av / 3)(sina, — cos a) and 
K' = — K, define the two valleys. The effective mass 
model of graphene is introduced by the ansatz 



tf(r) = VA [e lKr (j> (r + e 



-iKr 



(3) 
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where two pairs of envelope functions cj>± and 
correspond to the valleys K and K', are smooth on the 
scale a and A = (3\/3/2)a 2 is the area of the elementary 
cell. Arranging these four slowly varying amplitudes in 
a single vector 
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one arrives at the valley-symmetric Dirac Hamiltonian, 
Hq = —ihvcrV, where hv = 3ta/2 and <r = [a x ,a v ) is 
the vector of Pauli matrices in the sublattice space. 

The wave function on the honeycomb lattice (0) is ex- 
pressed as ^(r) = (u(r)|$) where 



(u(r)| 
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The representation of the T-matrix of the on-site pot 
tial impurity located at r is given by [l4| 

T = |u(r))V<«(r)| . 



;en- 
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The phase factors in Eq. ([5]) are responsible for a site 
classification, which depends on the orientation angle a. 
Using Eqs. ([5j H3) we represent the impurity T-matrix as 



T C = I (1 + (<T Z T Z 



JOc 
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where I = AV/2 is the scattering length, <j± = (a x ± 
ia v )/\/2, r± = (r x ± ir y )/V2, Ux,y,z and T X y Z are the 
Pauli matrices in the sublattice and valley space, respec- 
tively. The phases are 9± = ±a + 2Kr. Assuming an A 
site at r = 0, we encounter three possible values of the 
factor e 2iKr = e~ iKr = e 2 ™ c / 3 with c = 0,±1. Together 
with the sublattice index £ this yields six possible T- 
matrices for on-site impurities. In order to visualize this 
classification, let us introduce a color scheme by assigning 
six colors to lattice sites (three in each sublattice) corre- 
sponding to the phase factors e z9± = e ±la + 27Tlc / 3 , These 
colors define a superlattice (shown in Fig. [T]) with a pe- 
riod of three elementary unit cells and with six atoms per 
supercell. For a = 0, only three distinct phases remain, 
yielding the same three colors in both sublattices. 

In order to demonstrate the sensitivity of transport 
quantities to the colors of impurity sites, we compute 
the conductance of a rectangular graphene sample with 
two vacancies at positions ri and r?- The conductance is 
calculated numerically using a modification [15| of the re- 
cursive Green's functions algorithm [l6| that accounts for 
peculiarities stemming from the graphene lattice struc- 
ture, as well as analytically from the unfolded scattering 
formalism The numerical and analytical results are 
compared in Figs. [2J [3] 

We consider a rectangular graphene sample of the di- 
mensions LxW with metallic leads attached at x = and 
x = L. In order to model the pseudodiffusive transport 
regime, we assume W 3> L. In our analytic approach, 
the chemical potential in the leads is fixed infinitely far 
from the Dirac point, while the chemical potential inside 
the sample is set to zero. For two on-site impurities with 
scattering length £, the correction to the conductance as 
compared to the clean case is conveniently expressed as 
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where Gij — ( u ( r i)\S\ u (' r j)) are the matrix elements of 
the Green's function operator defined by 



/! — Hq + iO —a x r)S(x) 
—a x rjS(x — Li) fj, — Hq — iO 



t. 



(9) 



Here ij = i sinh(</>/2) is the counting field applied at the 
metal-graphene interfaces x = and x = L. Below we 
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study the case of two vacancies modeled by an infinite 
on-site potential. This corresponds to the limit I — > oo 
and allows us to simplify Eq. flSJ). 

The exact analytical solution to Eq. ^ with the spec- 
ified choice of the chemical potential was found in Ref. 
(loj . Using this solution together with the expressions 
§5§ in Eq. @ we express the conductance of a rectangu- 
lar graphene sample with two vacancies as G = Go + <5G, 
where Go = Ae 2 W/nhL is the conductance of the clean 
sample. The vacancies located at positions ri and 
inside the sample give rise to the correction 



4e 2 y 2 P11P22 ReO^) 



SG = -h 
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where y = y\ — 1/2 ■ Here we introduced the notation 

Pij = e l(e '-^ )/2 csc [ir(CiXi + (jXj+iyi-iVj)/2L\, (11) 

where 9i = aQ — Kr, is the T-matrix phase, which can 
take on one of six values 9±. The phase difference 9 — 
9\ — 62 encodes the dependence of conductance on the 
colors of vacancies and on the orientation angle a. 

For numerical simulations we have chosen the sam- 
ple length L « 50 a, the width W ~ 600 a, and con- 
sider two crystal orientations with a = (armchair) 
and a = 7r/6 (zigzag). The metallic leads at x < 
and x > L are defined by the large chemical potential, 
A*oo = 0.3 1 ~ lOHv/L (measured with respect to the 
Dirac point). Inside the sample (0 < x < L) the chem- 
ical potential is tuned to the close vicinity of the Dirac 
point, /l*o = 0.001 1 sa 0.033 hv/L. In the absence of va- 
cancies, the relative deviation of the conductance from 
the value Go = 4e 2 W/nhL is found to be less than 1%. 
While a single vacancy has only a negligible effect on the 
numerical results, a pair of vacancies leads to a finite 
correction, 5G, which we calculate numerically for dif- 
ferent positions of one vacancy while keeping the second 
vacancy fixed in the center of the sample. 

One set of data is shown in Fig. [3] for a = and for 
vacancies of different type: A- vacancy at ri = (L/2,0) 
and S-vacancy at r2 = (#2,2/2) where x-i ~ 2L/3. The 
dependence of 8G on all possible values of y = —yi 
on the lattice is shown with the data points. The con- 
ductance switches between three different smooth curves 
which correspond to a certain difference, 9 = 9\ — 02, of 
the T-matrix phase at the two vacated sites, 9 = and 
9 = ±2n/3, respectively. Similar data is plotted in Fig. 
[3K for the zigzag orientation of the crystal and for differ- 
ent positions of the _B-vacancy, X2 ~ 0.45L. Here three 
other phase differences, 9 = it and 9 = ±7r/3, appear. 
In Fig. [3p both vacancies are chosen to belong to the A 
sublattice (x2 ~ L/3, zigzag orientation). In this case we 
find again 9 = and 9 = ±2tt/3 in agreement with the 
classification of site colors presented in Fig. [TJ 

In Figs. [5J |3] we demonstrate a remarkable agreement 
between the numerical data obtained by the recursive 
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FIG. 2: (Color online) Conductance variation for an "arm- 
chair" (a = 0) sample with two vacancies (A and B). Chang- 
ing the distance y, the conductance jumps on the atomic scale 
between three different smooth curves corresponding to 8 = 
(green disks), 9 = 2n/3 (blue squares), and 9 = —2n/3 (red 
diamonds). The numerical data agrees well with Eq. (|10l) . as 
shown by the corresponding curves. 



Green's function technique and the analytical result (|10l) . 
The small deviations between analytical and numerical 
results are mainly due to the finite ratio W/L and the 
chosen finite chemical potential in the leads; a small de- 
tuning from the Dirac point within the sample, used in 
the simulations to avoid numerical instabilities, may also 
play a role. A good agreement with the theory is simi- 
larly obtained for other positions of the vacancies. 

Transport in the "short and wide" (W 3> L) sample 
can be described by the mean conductivity a = LG/W. 
For a low concentration of vacancies a virial expansion 
for a can be established. The lowest (second) order term 
of this series is readily found from Eq. pU|) : 



{i+i 4 Eh( CI -c 2 )( 



nA Cl nAc 2 +nB Cl nBc 2 ) 



4e2 
wh 



+ 2% a+ ^ (ci _ C2) n Aci n Bc2 \ +G(7i J i b )|, (12) 

where ua c and uq c with c = 0, ±1 are the impurity con- 
centrations of the corresponding sites of the graphene 
superlattice (Fig. []}. Performing numerical averaging 
over impurity positions in Eq. (|10p we find the values 
70 ~ 0.2653 and 7±2tt/3 ~ —0.1197, which are inde- 
pendent of the angle a. Both j x and 7 X are even 27T- 
periodic functions of x- The mean value of 7 X over the 
period is zero. In a sample with "armchair" orienta- 
tion (a = 0) relevant parameters are 70 ~ 0.1700 and 
7±2tt/3 ~ —0.0850. In a "zigzag" sample (a = 7r/6), we 
have 7±7r/3 ~ 0.0843 and 7^ ~ —0.1686. In each case, 
the sum of the three values is vanishingly small ~ 10~ 5 . 
If impurities are distributed uniformly among the lattice 
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FIG. 3: (Color online) Conductance variation for a "zigzag" 
(a = 7r/6) sample with two vacancies. Top panel: A and B va- 
cancies. On changing y the conductance jumps on the atomic 
scale between different smooth curves corresponding to 8 — n 
(green disks), 7r/3 (blue squares), and — 7r/3 (red diamonds). 
Bottom panel: Both vacancies in A sublattice. The phase 
difference 9 is either (green disks), 2-7r/3 (blue squares), or 
— 2-71-/3 (red diamonds). The numerical data agrees well with 
the result Eq. (flQ)l . as shown by the corresponding curves. 



sites, the conductivity is given by 

a= — (l + nn 2 L 4 + 0(n 3 L 6 )) 
Tin 



(13) 



where n is the impurity concentration and k w 0.0043 
exhibits tiny oscillations ~ 10 -5 as a function of a. 

However, the dependence of the conductivity on the 
crystal orientation is substantial if the distribution of 
adatoms among sites of different colors is not uni- 
form. Such a non-uniform distribution can originate from 
adatom correlations. For instance, the electron-mediated 
interaction of on-site impurities [l?], EH is inferred from 
the free energy F = —TJ2 U Tr\n(l — £(u\Q\u)), where Q 
is the Green's function © at Matsubara energy iu> n and 
(f> = 0. This gives rise to interaction oscillations yield- 
ing a color scheme similar to that depicted in Fig. [T] It 
is worth noting that Refs. [13] addressed the problem 
in an infinite system, whereas the effective interaction of 
adatoms with the leads [encoded in Eq. ©] may play an 
important role in a finite sample. 



In conclusion, we theoretically proposed and numer- 
ically confirmed an extended classification of impurity 
sites in the graphene honeycomb lattice for the case of 
strongly-bound adatoms or vacancies. The classification 
is illustrated in Fig. Q] by assigning colors to the lattice 
sites. The general analytical expression for the Dirac- 
point conductance of a graphene sample with two reso- 
nant on-site impurities is given as a function of impurity 
coordinates. The Dirac-point conductivity of graphene 
with a small number of randomly distributed adatoms 
is shown to be sensitive to the relative concentration of 
impurities at the sites belonging to different sublattices 
and having different colors. 
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